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Abstract: A method to calculate heavy flavor visible-level cross sections at Next-to-Lcading Order (NLO) in 
Quantum Chromodynamics (QCD), based on an interface of the FMNR program to Pythia, is described. It 
uses the NLO prediction at quark level provided by FMNR, with a statistical reduction procedure (REDSTAT) 
that allows a link to Pythia 6.2 to be made, from where the description of the full hadron fragmentation 
and decay chain is obtained. The method is applied to ep — » bbX — > D*fiX and ^j,^fi^X final states at 
HERA. Comparisons of the data and NLO cross sections at visible and 6-quark level were found to be consistent. 
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1 Introduction 

Heavy flavor (beauty and charm) production in ep 
collisions at HERA allows QCD to be probed and 
understood in detail. Several heavy flavor produc- 
tion channels have been studied by the HI [2] and 
ZEUS [1] experiments and compared to NLO QCD 
predictions based on the FMNR [3] massive fixed 
order calculations. These calculations are currently 
the only ones available at NLO which include a fully 
differential description of the parton- level final states, 
and therefore a non-trivial propagation of cuts at the 
visible level (i.e. on variables measured directly in the 
detector). In order to achieve this, the partons need 
to be fragmented and decayed to measurable hadron 
or lepton final states. For relatively simple cases, this 
can be implemented through a numerical or analytical 
parametrization of the parton-hadron level relation, 
obtained from separate MC simulations or direct 
measurements, which are then directly appended to 
the FMNR prediction. For final states with correlated 
cuts on several final particles, such as in [4, 5], this is 
not easily possible, although it has been attempted [6]. 
Nevertheless, NLO predictions for such final states are 
needed. 

One solution to this problem is the MC@NLO 
approach [7] - already implemented for the LHC and 
Tevatron - which is not yet available for HERA. In 
this approach NLO QCD predictions arc matched to 
parton showers in heavy flavor production, combining 
a Monte Carlo event generator (HERWIG [8]) with 



NLO calculations of rates for QCD processes. 

The alternative proposed here is an interface of the 
FMNR program to Pythia [9]. The application of 
the REDSTAT option [10] to the NLO description at 
quark level provided by FMNR, transforms FMNR 
into an effective Monte Carlo-like event generator. 
The parton-level events are then read back through 
the "Les Houches accord" [11] interface into Pythia 
6.2, where the description of the hadron fragmentation 
and full decay chain is implemented. 

No attempt is made to include additional parton 
showering. This will be achieved later by MC@NLO 
for HERA, when it will become available [12]. 



2 The FMNR program 

FMNR makes cross-section calculations at NLO in 
QCD for heavy quark production in ep and "/p colli- 
sions. It is designed to give predictions in the photo- 
production regime (photon virtuality < 1 GeV^). 
It contains point-like and hadronic photon coupling to 
the heavy quarks (see Fig. [T]) in the massive flxed-order 
approach. In this scheme, the quarks are produced 
dynamically (i.e. they are not included in the struc- 
ture functions) and the relevant scale is the mass of 
the heavy quarks (toq, Q — c,b), neglecting terms like 
[as In /^^/mq]" that are resummed in other approaches. 
This makes this approach particularly valid for pro- 
cesses in which the hard scale /i^ ~ TOq. 
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Figure 1: Point-like (left) and hadronic (right) photon 
coupling to the heavy quark in the leading order (LO) 
boson-gluon fusion process. 

In addition to quark-level predictions, the FMNR 
program provides a framework to fragment e.g. b 
quarks into B hadrons, and simulate the decay of 
these hadrons by interfacing them to appropriately 
chosen decay spectra. However, decays to complex 
final states, like D* + fi from the decay of the same B 
hadron with cuts on both particles, cannot be easily 
implemented in this scheme. 

A straightforward interface of the parton-level events 
produced by FMNR to a Monte Carlo-like fragmenta- 
tion and simulation chain is not practical since weights 
(either positive or negative) are assigned to the output 
events. These weights range over more than 8 orders 
of magnitude as is shown in Figure [21 FMNR also dif- 
fers from a true Monte Carlo generator because the 
events are not created randomly, but rather in a sys- 
tematic order, and thus the complete data set must 
be considered in order to get meaningful results. This 
makes such an approach extremely inefficient because 
high statistics needs to be generated in order to keep 
statistical fluctuations low. 

As seen in Figure [H in most cases the weights as- 
signed to the events are within four orders of magni- 
tude. However, there are also cases in which a pair of 
similar events are created with relatively large (but op- 
positely signed) weights. This is due to the fact that in 
the phase-space region where FMNR performs the cal- 
culations, coUinear and infrared divergences are present 
[13]. Physically, these paired events correspond to two 
distinct but similar processes: In the first quark- 
antiquark pair and a gluon are produced in the colli- 
sion, and the gluon is such that it remains near one of 
the heavy quarks. In the second case, the gluon not 
only remains near the heavy quark but also at some 
later time is reabsorbed (see Figure |3]). 
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Figure 2: An example of distributions of positive (con- 
tinuous line) and negative (dashed line) weights for 
events generated using FMNR. 
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Figure 3: Examples of processes with positive or neg- 
ative weight, and their corresponding contribution to 
the cross-section in the FMNR calculations. 



2 



For the cross-section calculation this means in 
practice that the process with e.g. the coUinear gluon 
has positive weight, but gives a divergent contribution, 
and the interference term between the LO process 
and the process with a virtual gluon has a negative 
weight and gives a divergent contribution as well. 
FMNR overcome the problem of the soft and collinear 
divergences by generating sequences of topologically 
correlated events in such a way that the cancellation 
is performed by the integration over that phase-space 
region. The output of FMNR is hence, weighted 
partonic events with a heavy quark-antiquark pair, 
and events with the pair plus an extra parton (a gluon 
or a light quark or anti-quark). 

However, problems do occasionally arise because 
the histograms used to plot the output are discretely 
binned, and it can happen that one of the two paired 
events will end up in one histogram bin, while the 
other event with only very slightly different topology 
ends up in the adjacent one. When this happens 
the two events can not cancel one another, and so 
the two adjacent bins are left with relatively large 
weighted events. It is then necessary to produce many 
more events with smaller weights for the histogram 
to smooth out. Although the contributions to the 
cross section are finite at this point, the values of 
the weights for the generated events span over many 
orders of magnitude (Fig.[2|). 

A solution to this problem is to find all paired events 
and suitably average them before they are output into 
the histogram. This is what the reduced statistics 
(REDSTAT) option does, as will be explained in de- 
tail in the following section. 



3 REDSTAT 

The group of subroutines REDSTAT is implemented 
as an extension to the FMNR program. It is designed 
to: 

• reduce the range of weights for the generated 
events, such that an interface to fragmentation 
and simulation packages like Pythia becomes re- 
alistic; 

• reduce the necessary statistics without losing NLO 
accuracy; 

• improve run-time efficiency. 

The method is as follows: REDSTAT monitors all 
the parton level events as they are generated by 



FMNR, and searches for events with weights above 
a given threshold. Two or more events with large 
but opposite sign weights and similar kinematics are 
then combined into a new single event by averaging 
the four-momentum of the partons and assigning the 
sum of the weights as the new weight. For events 
with weights below the threshold, REDSTAT makes a 
random decision to keep the event with a probability 
proportional to its weight (sampling approach). These 
events then acquire the weight of the threshold. 

The threshold is chosen such that the number of 
events with positive weights is almost equal to the 
number of events with negative weights above that 
threshold. Events are considered to have similar kine- 
matics when the difference in transverse momentum 
(pt), rapidity {Q and azimuthal angle (0) of the 
quark-antiquark pair are less than user cut values that 
should reflect the experimental detector resolution. 

After this procedure, as can be seen in Figure [U 
the range of weights is reduced to about two orders of 
magnitude, and the number of generated events is also 
reduced, therefore improving the run-time efficiency 
of subsequent applications using these events as input. 
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Figure 4: Distribution of weights for the generated 
FMNR events after REDSTAT. 

This method preserves the NLO accuracy for the 
relevant spectra (e.g. 6-quark px and angular distri- 
butions) and cross sections at parton level, as long as 
the chosen binning is not smaller than the assumed 
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resolution, as will be explained in the next section. 

Some caveats of this method are related to technical 
details of the REDSTAT implementation: the optimal 
threshold used to decide whether events will or will 
not be combined does not have a universal value. It 
must be chosen looking at the weight distributions, 
and these change depending on the generated process 
(charm or beauty) and according to the number of 
generated events (for example, it moves to the left 
when more statistics is generated). The default cut-off 
values of pt, C 4> are tuned to the ZEUS detector 
resolution and have to be appropriately retuned 
for other applications. Therefore, we recommend to 
contact the authors for advice before using REDSTAT. 

Finally, an optional output file for the interface to 
other packages is created. This is an ASCII file con- 
taining the information of the generated parton level 
events. For each event, this includes e.g. the assigned 
weight and the four-momenta of the two or three final 
state partons. 

3.1 Examples of REDSTAT results 

• Cross-section prediction at b-quark level 
for the process -fp b{b)X 

This was calculated using the original FMNR 
program and compared to the FMNR version 
which includes REDSTAT. The kinematic region 
was < I GeV^ and inelasticity 0.05 < y < 0.85. 

The parameters used in both calculations were: 

— the mass of the 6-quark nib — 4.75 GcV; 

— the renormalization and factorization scales, 
defined as /.j = -^/jti^ +-Pt6j where p^j^ is the 
average of the squared transverse momentum 
of the two emerging b quarks; 

— the structure functions CTEQ5M [14] for the 
proton and GRV-G-HO [15] for the photon. 

For the evaluation of the theoretical uncertainty 
(relevant for section O, the scales were varied by 
a factor 2, and the mass of the beauty quark was 
varied between 4.5 and 5 GeV. The variations 
were done simultaneously such that the spread 
was maximized. The uncertainty of the structure 
functions turned out to be small in comparison, 
and is neglected. 

Additional REDSTAT parameters were: 



— weight threshold: 10 ^; 

— maximal difference in transverse momentum 
between the partons of the events to be com- 
bined, ApT = 1.0 GeV; 

— similarly, maximal difference in rapidity: 
AC = 0.20; 

— maximal difference in the angle between the 
quarks: Acj) — 0.3 rad. 

These cuts were applied to the two heavy 
(anti)quarks only. The number of events gener- 
ated in each sample are displayed in Table [TJ 



Sample 


Number of events 


FMNR (original) 


31318299 


FMNR with REDSTAT 


3159440 


FMNR high statistics 


322391999 


FMNR fast REDSTAT 


3077217 



Table 1: Number of events generated per sample. 

The results shown in Table [2] demonstrate that 
the NLO accuracy of the original FMNR is 
preserved after the application of the REDSTAT 
option. 





FMNR (original) 


4.95 nb 


FMNR with REDSTAT 


4.94 nb 



Table 2: Comparison of NLO cross-section predictions 
using FMNR without/with the REDSTAT option, for 
rapidity of the 6-quark or antiquark |C| < 1. 

• Comparison of the transverse momentum 
distribution of the b(b)-quark 

As it is shown in Figure [5l the REDSTAT option 
yields a good description of the 5-quark px- The 
small differences in the shapes of the distributions 
are due to statistical fluctuations of such a typical 
FMNR run. 
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This is illustrated further in the FigurelHl in which 
in order to estimate the statistical fluctuations, 10 
samples of standard FMNR runs where generated 
using the same statistics changing only the seed of 
the random number generator. Obtaining a dis- 
tribution as the mean value of the 10 samples, the 
error of which has a variance V{x) = a'^/N. This 
mean distribution (called standard in the Figure) 
is compared to a run using 10 times more statistics 
and to the version using REDSTAT. After RED- 
STAT, the fluctuations arc expected to be smaller 
than in the original FMNR, since the probabil- 
ity to split a correlated pair of events at the bin 
boundary is reduced. 




2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 

Pt (GeV) 



Figure 5: Comparison of the 6(6)-quark transverse mo- 
mentum distributions using the original FMNR (con- 
tinuous line) and FMNR with REDSTAT (dashed 
line). 



Variation of detector resolution parameters 

A comparison between a run using default values 
tuned for ZEUS {Apr = 1.0 GeV, AC = 0.2 and 
A(/) — 0.3 rad) and a run with wider cut values 
(ApT = 1-5 GeV, AC = 0.38 and A0 = 0.45 
rad) which therefore allows a faster performance, 
is shown in Figure [71 The fluctuations result from 
the differences in the event combination proce- 
dure, since the same original FMNR events are 
generated in each case. 
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Figure 6: 5(6) transverse momentum distributions 
of the original FMNR with 10 times more statistics 
(squares), a mean of 10 FMNR samples (full circles) 
where the error bars are the statistical errors properly 
accounted for the different event weights, and FMNR 
using REDSTAT (stars). 
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Figure 7: Comparison of the 6(6)-quark transverse mo- 
mentum distributions using REDSTAT with default 
parameters [1 GeV, 0.2, 0.3 rad] and REDSTAT with 
wider cut values [1.5 GeV, 0.38, 0.45] (fast). 
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4 The FMNR(g)PYTHiA Interface 

The implementation of the interfaee of the FMNR 
parton-level predictions to the fragmentation and 
decay chain from Pythia/Jetset is a two-step pro- 
cess. The first step of this FMNR®Pythia interface 
consists in the apphcation of the REDSTAT extension 
to the FMNR program (described in Section RED- 
STAT is used to transform FMNR into an effective 
Monte Carlo-hke parton-level event generator, through 
the combination of events with similar kinematics and 
a sampling approach. 

The events obtained this way are written to an 
output file. In the second step, these events are read 
back into the Pythia 6.2 program through the "Les 
Houches accord" user interface. 

The interface routine UPEVNT then sequentially 
reads one event from the input ASCII file and fills it 
into the "Les Houches" interface variables to be used 
as input by Pythia. At this stage, a "reasonable" 
(i.e. physically possible) colour fiow is assigned to each 
FMNR parton level process. FMNR does not provide 
this information, which is needed in the case of string 
fragmentation. The difference due to different possible 
colour fiow assignments should be included in the 
systematic error if this is critical for the application. 

Since the first comparisons with data were made 
to measured beauty production cross sections [], the 
following sections are written for the specific case of 
beauty production. All arguments hold in an analog 
way for charm production. 

The initial state partons are allowed to have an 
intrinsic kx (typically ~ 300 MeV) as implemented in 
Pythia. This has a negligible effect on the resulting 
cross sections (~ 1%). 

Parton showering is not allowed in order to avoid 
double counting of higher order contributions. 

Fragmentation of 6-quarks close to production 
threshold turned out to be a non trivial issue. Since 
the details of the threshold treatment were found to 
be much more important than the choice of a particu- 
lar fragmentation function, the Peterson formula with 
e = 0.0035 is used for convenience. Three approaches 
have been considered: 

• Independent fragmentation in the PYTHIA model, 
within FMNR(g)PYTHiA. 

This is used because FMNR does not provide color 



connections on an event-to-event basis, and color 
connections are not required in this model. 

• Fragmentation in the Lund string model, within 

FMNR(g)PYTHIA. 

For this, reasonable color connections have to be 
associated to each FMNR event. 

• Independent fragmentation scheme as provided in 
the context of the original FMNR. 

This implies setting the B-hadron momentum 
equal to the 6-quark momentum before reducing it 
according to the Peterson formula. This neglects 
threshold corrections due to the need of simulta- 
neous conservation of energy and momentum. 

In the first two cases, the fragmentation, decay tables 
and kinematics as implemented in Pythia 6.2/Jetset 
are used to obtain a hadron-level event. Therefore, 
non-dominant arbitrarily complicated decays, such as 
B ~> D*D followed by a 13 ^ jiX, or muons through 
intermediate J/"^ or r states, are automatically 
included. The branching ratios were empirically 
corrected at analysis level to correspond to those 
obtained from the Particle Data Group (PDG) [16]. 
Unless otherwise stated, the 2"'' approach is used 
for the central prediction, while the first enters the 
systematic uncertainty. 

The third case could be used only for relatively sim- 
ple final states (e.g. independent production of a muon 
from each heavy quark), and the resulting differences 
with respect to the first two cases were extrapolated to 
the more complicated ones (e.g. correlated production 
of two muons from the same b quark), and used as a 
systematic check. 

5 Examples of the use of the 

FMNR(g)PYTHiA interface. 

Double tagging techniques like fi^n" or D*iJ. [4-6] 
benefit from the increasing amount of data provided by 
HERA. At the ZEUS experiment, a beauty enriched 
data sample is obtained in these channels taking 
advantage of the significantly reduced background. 
This allows softer kinematic cuts which, in combi- 
nation with the wide rapidity coverage of the ZEUS 
muon detectors, results in an enhanced sensitivity 
to B hadrons produced at low transverse momenta, 
where the bulk of the b cross section is concentrated. 
Under these conditions a direct measurement of 
the total bb cross section becomes possible without 
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the use of large model-dependent extrapolation factors. 

Tagging muons or D* from different b quarks allows 
bb correlations to be explicitly measured, while tagging 
two muons or a muon and D* from the same b quark 
yields measurements which arc almost insensitive to 
such correlations. However, the complexity of the 
decays of these channels and the different cuts applied 
to the final state particles makes a calculation of the 
NLO prediction using the original FMNR approach 
very difficult. 

In this section our new method will be applied to 
obtain NLO beauty cross-section predictions for the 
final state channels D*fi and ^~^fi~ . 

5.1 The ep hbX D*^iX' channel 

In the case of D* + /i [4] an unlike-sign combination 
of both particles is observed in the same detector 
hemisphere if they come from the same b quark, 
(mainly via B*^ D* ^v^) yielding a pure b sample. If 
they originate from different b quarks then the /i and 
D* are emitted into different hemispheres. Both like 
and unlike sign charge combinations are allowed. In 
the unlike sign case only, the signature is similar to 
the one given by charm from different c quarks. These 
differences can be used to extract the beauty signal. 
More details are given in [17]. 

The cross section for inclusive beauty production 
ep bbX D*fiX' in the kinematic range > 1.9 
GeV, -1.5 < r;^' < 1.5, > 1.4 GeV and -1.75 < 
77^ < 1.3 was measured to be: 

cr„,^ = 160 ±37(stat.)t^?(syst.) pb. (1) 

This is larger than, but still compatible with, the 
corresponding FMNR(8)Pythia NLO prediction of 

al^o = 67t?;(NLO)l^^(frag. © br.) pb. (2) 

where the first error refers to the uncertainties of the 
FMNR parton level calculation, and the second error 
refers to the uncertainties related to fragmentation 
and decay. 

Strictly speaking, the FMNR predictions are only 
valid for the photoproduction regime. Here, the 
Weizsacker- Williams (WW) approximation with an 
effective Ql,^^ < 25GeV^ cutoff [18-20] was used to 
include the ^ 15% DIS contribution for a combined 
cross section. This NLO QCD prediction is also listed 
in Table [3 



A cross section for the same kinematic range, but 
adding a photoproduction requirement (Q^ < 1 GeV^ , 
0.05 < y < 0.085) was also obtained. The result, as 
well as the corresponding FMNR(X)Pythia prediction, 
are shown in Table [H As in the inclusive case, the 
prediction underestimates the measured cross section, 
but is compatible with the measurement within the 
large errors. 

Finally, these measured visible-level cross sections 
were extrapolated to &-quark level using Pythia. The 
NLO QCD prediction can then be obtained at par- 
ton level directly from the original FMNR calcula- 
tion. From the comparison of the ratios at visible 
and 6-quark level in Table [3 one can conclude that 
the Pythia extrapolation, which was the only way to 
compare the cross sections before Eq. ([2|) could be cal- 
culated [17], was reliable. 

5.2 Comparison ZEUS - HI 

The HI Collaboration has measured a cross sec- 
tion similar to the ZEUS photoproduction cross 
section, in a slightly different kinematic region [6]: 
p^' > 1.5 GeV, -1.5 < jf < 1.5, > 2.0 
GeV, -1.735 < rf < 1.735, < 1 GeV^ and 
0.05 < y < 0.75. Its value is shown in Table HI Using 
the FMNR®Pythia interface, the ZEUS cross section 
can be extrapolated to the HI kinematic region [4]. 
A direct comparison of the two results is shown in 
Table m Reasonable agreement is found. 



Table 4: Comparison of HI and ZEUS visible cross sec- 
tion for ep — > bbX — > D*^X'. The ZEUS cross section 
has been extrapolated to the HI kinematic range. 

cross-section HI ZEUS (cxtrap.) 

7P 206 ± 53 ± 35 pb 135 ± 3313^ pb 



The corresponding NLO prediction from 
FMNR(g)PYTHiA is 

<yv^s.Hl = 6lll^(NL0)+f (frag. © br.) pb. (3) 

The data to NLO ratio is again consistent with the 
ones quoted in Table [31 

The NLO prediction is somewhat larger than the one 
evaluated in [6] due to the inclusion of the hadron- 
like photon contribution, the inclusion of secondary- 
muon branching fractions for D* and /i from the same 
b quark (e.g. B D*D D* jiX), which are difficult 
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Table 3: Comparison of measured and predicted cross sections. For the measured, the first error is statistical 
and the second systematic. The NLO predictions at visible level were obtained using FMNR(g)PYTHiA, and the 
measurement at &-quark level was obtained extrapolating with plain Pythia. 



cross-section 




measured 


NLO QCD 


ratio 


Visible 


inclusive 

IP 


160 ± 37j:^? pb 
115 ± 29+27 pb 


67t?l pb 
54^i« pb 


2 4+0.9 

1 +0.8 

^■-•--i.o 


b level 


IP 


11.9 ±2.9^3:^ nb 




2.ot;:« 



to handle outside the FMNR(giPYTHiA framework, and 
a detailed simulation of the kinematics of the b B ^ 
D* chain rather than direct coUinear fragmentation of 
b quarks into D* mesons. This reduces the discrepancy 
claimed in [6] and is an example for the importance of 
a detailed treatment of complicated final states. 



5.3 The ep ^ bbX 



X' channel 



The signature of the dimuon channel is very similar to 
the D* n channel and detailed information can be found 
in [5]. Again, this is a very complicated final state, for 
which the cross section measurements are being final- 
ized by ZEUS. The purpose of this section is to provide 
cross section predictions to which the ZEUS results can 
be compared. 

5.3.1 Visible beauty cross-sections prediction for 

ep ^ bbX ^l+^l-X' 

The measured visible cross-section compared to the 
NLO prediction from FMNR®Pythia is shown in 
Table [H The extrapolation to &-quark level was done 
using Pythia and is compared to the NLO prediction, 
obtained by adding up the (original) FMNR and 
HVQDIS [21] predictions, for the photoproduction 
and DIS regions respectively. 

Here, as in the case of the D* ji channel, the cross- 
section comparisons at visible and 6-quark level are 
consistent. 



ZEUS 



> 

a> 
O 

Q. 



Q. 



10 



10 



S" 


45 






i 


40 




35 


"B 


30 


■o 




25 




20 




15 




10 




5 








• ZEUS (prel.) 96-00 

(PYTHIA+RAPGAP) x 1 .95 

Pq NLO QCD (FMNR*PYTHIA) 



ZEUS 



9 10 

PT^'(GeV) 



• ZEUS (prel.) 96-00 

(PYTHIA+RAPGAP) x 1 .95 

'EB NLO QCD (FMNR*PYTHIA) 



-1.5 -1 -0.5 0.5 1 1.5 



2 2.5 



5.3.2 Differential cross sections 

Differential cross-sections were also obtained and 
shown in Figure [Sj One can observe the general trend 
of NLO to lie below the data, but consistent between 
errors. The shape of the distributions is well repro- 
duced by the NLO prediction. 

For the correlations between the two 6-quarks, the 
reconstructed dimuon mass range was restricted to 



Figure 8: Differential cross-sections da/dpr (left) 
and da/drj (right) for muons. The data points 
are compared to the scaled LO prediction from 
Pythia + Rapgap, and to the NLO prediction from 

FMNR(g)PYTHIA. 
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Table 5: Comparison of measured and predicted dimuon cross-sections. For the measured, the first error is 
statistical and the second systematic. 





cross-section 


measured (prel.) 


NLO QCD 


ratio 


Visible 


Total 


63 ± 7tTs pb 


30^^ pb 


9 1+0.8 


b level 


Total 


16.1 ± 1.8t^;? nb 


6.8if ? nb 


9 q+l-O 



m^'^ > 3.25 GeV. This additional cut reduces the prob- 
ability that both quarks come from the same B-hadron. 
The corresponding differential cross-section is shown in 
Figure [SI The distribution is reasonably described by 
the FMNR(X)Pythia NLO prediction within the large 
errors. The leading order contribution alone, which 
is also shown, reproduces the measured less well. Al- 
though the difference in shape is not dramatic, this 
could confirm the importance of the contributions from 
higher order processes. 
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Figure 9: Differential cross-section (ia/A0 for muons 
from different &-quarks 



5.4 Other applications 

In principle, the FMNR(X)Pythia interface can be 
used to generate events which can be propagated 
through the full simulation and reconstruction chain 
of any experiment measuring ep or heavy flavor 
final states (if the chain allows the handling of negative 
weights). This opens up additional applications, such 
as cross sections for hadron-level heavy flavor jets, 
with hadron to parton-level corrections fully consistent 
with the NLO approach. So far. LO+PS Monte Carlos 
needed to be used for this purpose. 



Examples for potential applications in the near fu- 
ture include 

• evaluation of parton to hadron-level corrections for 
measurements of the charm fragmentation func- 
tion at HERA [22, 23]; 

• evaluation of NLO cross sections for measured /i/i- 
jet-jet and /ie-jet-jet final states from both charm 
and beauty [24, 25], for which no NLO predictions 
exist so far; 

• reevaluation and verification of NLO predictions 
for essentially all earlier HERA heavy flavor re- 
sults [1, 2] at visible level. 

These applications will not obliviate the need for cal- 
culations using the MCQNLO approach, which will be 
able to handle final state parton showering in addition. 
However, they have the virtue of being available imme- 
diately. 

6 Conclusions 

Cross sections for heavy quark production in ep 
collisions at HERA are being measured with high 
precision. This includes channels with complex final 
states for which NLO QCD predictions are not easily 
obtained from simple extensions to parton level cal- 
culations. The FMNR(g)PYTHiA interface allows such 
predictions to be obtained for the photoproduction 
regime. 

In the ep — > hhX D* jiX' channel, both at 
visible and quark level, the measured data at ZEUS 
exceed the FMNR®Pythia NLO prediction, but are 
compatible within the errors. It was also found that 
the HI and ZEUS visible cross sections are consistent. 

The almost constant data to NLO ratio obtained 
from the cross sections shows that the extrapolation 
from visible to quark level, or vice versa, is meaningful 
and reliable. 
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Inclusive and differential cross sections in the ep 
bbX — > fi'^ii~X' channel are evaluated, which can be 
compared to ZEUS measurements, currently being fi- 
nalized. Other potential applications to both charm 
and beauty measurements in ep and 7p interactions 
have been outlined. 
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